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ABSTRACT 

FU Orionis systems are young stars undergoing outbursts of disc accretion and where the 
optical spectrum contains lines associated with both the disc photosphere and a wind compo- 
nent. Previous observations of the prototype FU Orionis have suggested that the wind lines 
and the photospher ic lines are modulated with periods of 14.54 and 3.54 days respectively 
(Herbig et al. 2003). We have re-observed the system at higher spectral resolution, by mon- 
itoring variations of optical line profiles over 21 nights in 2007 and have found periods of 
13.48 and 3.6 days in the wind and disc components, consistent with the above: this implies 
variability mechanisms that are stable over at least a decade. In addition we have found: i) 
that the variations in the photospheric absorption lines are confined to the blue wing of the 
line (centred on velocity ~ — 9 km s _1 ): we tentatively ascribe this to an orbiting hotspot 
in the disc which is obscured by a disc warp during its receding p hase, ii) The wind p eriod 
is manifested not only in blue-shifted Ha absorption (as found by Herbig et all (2003)) but 
also in red-shifted emission of Ha and H/3, as well in blue-shifted absorption of Na I D, Li I 
and Fe II A5018. iii) We find that the periodic modulation of blue-shifted Ha absorption (at a 
velocity of around —100 km s _1 ) is phase lagged with respect to variations in the other lines 
by around 1.8 days. This is consistent with a picture in which variations at the wind base first 
affect chromospheric emission and then low velocity blue-shifted absorption, followed - after 
a lag equal to the propagation time of disturbances across the wind's acceleration region - by 
a response in high velocity blue-shifted absorption. Such arguments constrain the size of the 
acceleration region to be ~ 10 12 cm. We discuss possible mechanisms for periodic variations 
within the innermost 0.1 AU of the disc, including the possibility that these variations indicate 
the presence of an embedded hot Jupiter in this object. 

Key words: accretion, accretion discs - line: profiles - stars: individual: FU Orionis - stars: 
mass-loss - stars: pre-main-sequence - stars: winds, outflows 



1 INTRODUCTION 

Although FU Orionis systems only represent a small class of 
young stellar objects, the benefits of understanding the physical 
processes present in these systems have implications for the whole 
of not only star, but also planet formation. FU Orionis systems 
are unique objects which provide the only opportunity to directly 
observe an accretion disc around a young star. Originally identi- 
fied as a sub-set of p re-m ain-sequence sta r s ent itled 'FUors' by 
Ambartsumianl dl97 ll) and lAmbartsumvarj J 197 lh . they are char- 
acterised by a large increase in optical brightness of ~ 4 mag- 
nitudes or more over a period of 1—10 y ears, with a slow de- 
cline of 50-100 years jHerbi d 1 1 96d 1 1 9771) . Other defining prop- 
erties of the class include spectral types of late F to G in the 
optical, with infrared spectra showing K-M giant-supergiant type 
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features dMould et alj|l978l: lHartmann & Kenvor]|l985l . ll987iK 
iKenvon et alJl988l : lHartmann & Calvetl 19951) . FU Orionis systems 
are associated with young T-Tauri stars, ev ident from the strong 
Li I A6707 absorption observed by iHerbid dl96fj) . a characteris- 
tic of young stars. The discovery of the second class member, 
VI 057 Cyg, f or which a pr e-outburst spectrum resembled that of 
a T-Tauri star dHerbid l 1958) supported this association with young 
stars. All FU Orionis objects are associated with star forming re- 
gions, as is evident from the ring shaped reflection nebulae first 
noted by iGoodrichl dl987l) . from the presence of infrared excess 
due to circumstell ar dust, and b y their general location in heavily 
extincted regions lHerbisj [l977h . More recent developments have 
suggest ed these objects tend to have binary companions, (see ta- 
ble 1 in I Vittone & Erricol {2005)), although the binary fraction is 
not necessarily higher than in other young st ars. A comp a nion c an- 
didate to FU Orionis itself was detected by Wang et al. (20 041) at 
a sepa r ation of 0.5 arcsecs. See Hartmann et al. I dl993h : lKenvonl 
dl995h : lHartmann & KenvorJ dl99fj|) : lHartmannl dl998l) for reviews. 
FU Orionis outbursts are thought to be produced by long pe- 
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riods of rapid accretion dLarsonlll983l ; E in & Papaloizou 1985), 
when the accretion rate is so great, ~ 10 _4 MQyr _1 dKenvon et alj 
1988), the inner disc heats up through viscous dissipation and the 
mid plane temperature exceeds that of the photosphere. Under these 
circumstances, the disc produces an absorption line spectrum with 
lower surface gravity than the central s tar, and the star's luminosity 
is swamped b y factors of 100-1000 lHartmann & Kenvon 1996). 
Conse q uently. lHartmann & Kenvonl 1 1985 , 1987a ); Kenvon et alj 
dl988h : IWeltv et alJdl992l) : |Popham et al.ldl996h : IZhu et al.ld2009h 
showed the spectrum of FU Orionis is well fitted by the pre- 
dictions of a broad band steady state accretion disc model. The 
periods of rapid accretion are thought to be triggered by ther- 
mal instabilities in the disc associated with partial ionisation of 
disc material, which ca uses the release of damme d up material 
on to the central star ( Hartmann & Kenv onl I1985D. The trigger 
mecha n isms are discussed by Pringle ( 1976); Hartmann & Kenvon 
1985 1) : iLin & Papaloizoul dl985h : Idarke et alj dl990h : iBell etail 



lysoj; Lin & rapaioizou ayaD), nance et ai. iiwuj; eeii ei ai. 
1995): Idarke & SveJ dl996l) : lLodato & Clarke! d2004l) : lzhu et alj 



( 2009), and references therein. The frequency of FU Orionis events 



observed within the stellar neighbourhood has led to the belief 
that many low mass stars may undergo mu ltiple FU Orionis events 
lHerbigj |l977l :l Hartmann & Kenvonl Il985l) . It is currently thought 
that up to 5 — 10 per cent of the matter accreted on to young T-Tauri 
stars occurs thro ugh repetitive periods of rapid accretion in FU Ori - 
onis outbursts ((Hartmann & Kenvonlll987al ; lHartmann et alj|l993l) . 
This highlights the importance of understanding this episodic ac- 
cretion process in the context of stellar formation. 

The high accretion rate in FU Orionis objects is accompa- 
nied by high mass-loss rates with most objects associated with 
molecular outflow s, optical jets and HH objects ( iReipurthl 1 1 990l : 
lEvans et "alii 19941) . The mass-loss rate i s around ~ 10~ 5 M yr" 1 
in FU Orionis itself dCroswell et all 1 19871 ; ICalvet et alj 1 19931 ; 
lHartmann & Calvetlll995l) . The strong bipolar outflow has a termi- 
nal wind velocity of ~ 300 km s" 1 , which has been detected pre- 
vi ously in the deep, broad blueshifted co mponents of Ha and Na I 
D jHerbid[l977l ; lBastian & Mundtll 19851) . It is thought the outflow 
arises from the surface of a Keplerian disc, ( lHartmann & CalveM 
1 19951), with the possibil ity of a disc 'chromosphere' first suggested 
bv lCroswell et alj dl987h to explain t he detection of Ha em ission, 
and modelled in the line profiles by iD'Angelo et al. | fcOOd) . The 
mechanism producing the rapid acceleration of material into these 
powerful outflows in circumstellar discs is not entirely clear. Possi- 
ble models involve the accretion flow spinning up the star close to 
break up and launching magneto-cent rifugal winds a long opened 
up field lines anchored in the star, dShu et al"] [l994). or conical 
winds driven by the pressure gradient of the azimuthal magnetic 
field component an d launched along opened up stellar field lines, 
( Koni gl et al1l201 ]]). However magnetic fields anchored in the ro- 
tating disc itself coul d also centrifugally accele rate material out- 
wards as described by Blandford & Payne ( 1982). The detection of 
a ~ 1 kG magnet ic field directed tow ards the observer at 0.05 AU 
in FU Orionis bv lDonati et alj d2005h supports the involvement of 
magnetic fields, yet their topology and role in driving mass-loss 
remains unclear. 

In addition to the variability on long time-scales, it appears 
FU Ori also exhibits unexplained variability o n shorter time- scales. 
Variability in Ha emission was first noted bv lHerbid l fl966T) over 
a time-scale of one y ear, w hich was confirmed in absorption and 
emission bv lHerbid (fl977h. This variability was observed to ex- 
ist on_a_jnoirthly_jime2Scale in observations of Ha and Nal D 
by iBas tian & Mundt ( 1985), who concluded the mass outflow ac- 
tivity from FU Orionis objects was a long lasting active state. 



This was confirmed bv lCroswell et all dl987l) who went on to sug- 
gest the maximum size of the region producing the variation ob- 
served in the absorption profiles of Ha at ~ 300 km s _1 was 
rough ly twice the optica l phot ospheric radius, or ~ 3 x 10 12 
cm. Hartmann & Kenvonl Jl985h noted variability in the blue spec- 
tral region and in the Mgl A5183 line and attributed the cause to 
variation in the mass ejection, concluding the wind arises from 
the optical disc surface in FU Orionis. Random, small-amplitude 
photometric variatio ns on time-scales of 1 day were observed by 
IKenvon et al] d2000l) . and the rapid variations were attributed with 
the flickering of the inner accretion disc in FU Orionis. 

The presence of periodicity was first evident in the bright- 
ness fluctuations a ft er the outburst, where Fourier analysis by 
IChochol & Tremkd dl980l) detected periods of 1645 and 3760 
days from observations between 1936-1976. Smaller cyclic vari- 
ati ons were then detected f rom photometry taken in 1984-1985 
bv lKolotilov & Petrovl ( fl985h . which showed a brightness period- 
icity in the V-band of 18.35 days for two months in 1984 and 
interestingly, a periodicity of approximately half, 9.19 days for 
the remai nder of the obser vations. The shorter period was con- 
firmed by Ibragimov ( 1993]), who in addition detected a period of 
9.2 days in the F-band li ghtcurve from 1987 August - 1988 Jan- 
uarv. lKenvon et al] d 19881) highlighted variability in FU Orionis be- 
low 4000A, which maybe attributed with the periodicity found by 
Kolot ilov & Petrovl fl985), or could be due to fluctuations in air- 
mass. Further confirmations of periodicities were found in the high 
resolution spectroscopy (13 km s -1 ) of 20 nights observed from 
1997-1999 analysed bv lHerbig etail d2003l) . who not only found 
variation in Ha on daily time-scales, but also confirmed a period 
of 13 to 18 days, which peaked at 14.847 da ys in the equ i valent 
width of Ha between -110 and -270 km s- MErrico etail d2003l) 
attributed this period to a rotationally modulated wind structure, 
envisaging a non-axisymme tric magnetic s tructu re co-rotating with 
the star. Note that whereas He rbig et alj d2003l) found a periodic 
signal only in the absorption component of Ha, Ivittone & Erricol 
(2005) detected periodicity in the Ha emission component of 6.70 
days in 1998-2000, but found no correlation and no periodicity 
present in the absorption component of Ha. 

It is not only the strong blueshifted absorption l ines affected 
by ma ss-loss that show periodicity in FU Orionis. iHerbig et alj 
d2003l) analysed spectral regions of weak photospheric lines, with- 
out strong blueshifted components (5540-5640 and 6320-6440 A) 
through cross correlation with a GO lb reference star (/3 Aqr) over 
20 nights between 1997 to 1999. This analysis demonstrated that 
the radial velocity of the peak of the cross correlation function 
undergoes periodic variations with a period of 3.542 days. It 
was suggested that this periodicity could be attributed to a non- 
axisymmetric hotspot, dHerbig et alj|2003l ; lErrico et alj|2003h . In 
the absence of any mechanism to sustain this hotspot against shear, 
a hotspot of radial extent 8R would be sheared out over R/8R or- 
bits. Thus it is of interest to discover whether this periodicity is 
sustained over many years. Long-lived periodic behaviour would 
be consistent, for example, with non-axisymmetric disc structure 
induce d by a stable magnetic co nfiguration or even an embedded 
planet dClarke & Armitagdl2003l) . 

In this paper we present observations taken over 21 nights 



1 Variability also appears to be present in the 6170 A spectral region of FU 
Orionis, which contains weak absorption lines from t he photosphere of the 
disc, however no time series analysis was conducted I Hart mann & Kenvonl 
11987 J) . 
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from the SOPHIE echelle spectrograph spanning a time frame of 
38 nights in 2007. The dense time coverage and the high resolution 
of 4 km s" 1 has allowed the periodicities in FU Orionis to be stud- 
ied in great detail. We observe the stability of the periodicitie s due 
to non-axisymmetric structure detected by Herbig et aL ( 20031) over 
a decade and highlight the necessity of a mechanism which main- 
tains such a structure against shear. In section [2] the details of the 
observations and data reduction techniques are described, with the 
variability in the line profiles described in sections [3] [4] and [5] The 
periodicity in the cross correlations of spectral regions containing 
weak photospheric lines is examined in section [6] Possible mech- 
anisms to sustain the periodicities over decades are discussed in 
section|7] 



2 OBSERVATIONS AND DATA REDUCTION 

Optical spectra of FU Orionis were obtained from the SOPHIE 
(Spectrograph pout 1' Observation des Penomenes des Interieurs 
stellaires et des Exoplanetes) cross-dispersed echelle spectrograph 
mounted on the cassigrain focus of the 1.93-m telescope at the Ob- 
servatoire de Haute-Provence (OHP, France), covering the wave- 
length range 3872-6943 A. The observations were obtained using 
the high resolution mode (R= 75000), over 21 nights, providing 
an excellent opportunity to study the optical line variability of this 
object in detail. A reference spectrum of f3 Aqr was also taken 
from SOPHIE in the same mode and used as the template in the 
cross-correlations. Th is GOIb supergiant is used a s the standard as 
it has been shown by lHartmann & Kenvorj J 1985f) that G spectral 
types provide a better signal-to-noise correlation than K spectral 
type templates, and it has the same effective temperature as the re- 
gion of the d isc of FU Orionis vis ible at optical wavelengths, as 
described by ICroswell et al.l dl987l) . The data were reduced auto- 
matically by the pipeline, which includes bias subtraction, optimal 
order extraction, cosmic ray removal, corrections for flat fielding 
and wavelength calibration. The 39 echelle orders were divided by 
the blaze function, reconnected and barycentrically corrected, re- 
sulting in a re-sampled spectrum with a constant wavelength step 
of 0.5 km s _1 , smaller than the original spectral sampling. Details 
of the observations are shown in Table Q] 

The data were continuum subtracted by clipping out the spec- 
tral features and performing median and boxcar filtering over a 
scale length of 20A, to obtain a pure continuum spectrum, which 
was then subtracted fr om the featured spect rum and normalized 
to unit continuum. See Batta glia et al] (2008) for more details on 
this method. This was performed over the wavelength range 4100- 
6865 A to ensure variations in the continuum were not included 
in the resultant line profiles and to minimize the propagation of er- 
rors dues to edge effects. Fig. [T] shows the initial normalisation of 
the JD=2,454, 106.309 observation over this region. Variance in the 
continuum was minimized by averaging all 21 nights of the nor- 
malised spectra and subtracting this from each spectra individually 
to leave the residual spectra. These residuals were then continuum 
subtracted again and the result was re-combined with the initial av- 
erage spectrum to ensure fluctuations in the continuum do not affect 
the line profile variability. 

Spurious emission was removed by flagging the points which 
were above three sigma from the mean, both in each spectra indi- 
vidually, and across simultaneous points in all 21 observations, and 
replacing with the result of a median filter across a region of ±0.1 
A around the bad value. This ensured large variations in emission 
over the spectral series were removed, reducing the potential con- 



Table 1. Observations of FU Orionis and /3 Aqr. 



Object 


Date 


Sequence 


Exposure 


Barycentric JD 






ID 


Time (s) 


-2,454,000. 


FU Orionis 


03-01-2007 


220020 


1800.0 


104.31945010 




04-01-2007 


229700 


2000.0 


105.33395783 




05-01-2007 


235720 


2500.0 


106.30874292 




06-01-2007 


244940 


2200.0 


107.30133838 






255950 


777.0 


107.32028452 




08-01-2007 


258300 


2000.0 


109.31690541 




12-01-2007 


274990 


2000.0 


113.31687568 




13-01-2007 


275710 


2000.0 


114.31604553 




14-01-2007 


276540 


1800.0 


115.28749519 






276550 


1800.0 


115.31044553 




15-01-2007 


277830 


2000.0 


116.37177376 




17-01-2007 


279210 


2000.0 


118.31951493 




18-01-2007 


281980 


2000.0 


119.42023451 




in A i n \~7 

ly-0 1-2007 


282860 


2000.0 


120.3 1694694 




20-01-2007 


284260 


2000.0 


121.33335079 




26-01-2007 


294910 


2000.0 


127.37686294 




27-01-2007 


295440 


2000.0 


128.30713983 




28-01-2007 


296290 


2000.0 


129.30472606 




29-01-2007 


297200 


2000.0 


130.29181404 




30-01-2007 


298190 


2000.0 


131.30359823 




31-01-2007 


299070 


2000.0 


132.30679379 




05-02-2007 


303430 


1800.0 


137.33476702 




10-02-2007 


306160 


1800.0 


142.33284243 


/3Aqr 


08-10-2007 


465830 


300.0 


382.34761084 



tribution from spurious lines. Telluric lines were identified through 
visual comparison between all the unbatycentrically corrected spec- 
tra stacked with the median of all spectra. Aligned, narrow lines 
which persisted in the median spectrum were clipped and removed 
from further analysis. 

The continuum subtracted spectra were cross-correlated with a 
template GOIb reference star, namely f3 Aqr taken from the sa me in- 
strument, via the methods described in iTonrv & Davisl dl979h . The 
template was shifted to zero velocity, as it has an radial velocity 
measured to be 7.11 km s _1 , when cross correlated with a syn- 
thetic spectrum of a GO star with T e g = 6000 K, log(g)= 3.0. 
This specific model was chosen as it provided the best match to 
the observed line strengths in the /3 Aqr spectrum. The radial ve- 
locity was derived by fitting a Gaussian to the peak of the symmet- 
ric cross-correlation function. This was necessary as previo us mea- 
surem ents were made with far less accurate resolution by IWilsonl 
dl953h . The observed template was continuum subtracted in the 
same manner as described above. Furthermore, the synthetic stellar 
spectrum was smoothed using a Gaussian with full width half maxi- 
mum (FWHM), comparable to the corresponding unblended, unsat- 
urated absorption lines in the relevant regions of the /3Aqr spectrum. 
The resulting spectrum was used as an addition reference standard 
in cross correlations in order to confirm the zero point of the spec- 
tral lines in the /3 Aqr spectrum. Using this method over a spectral 
interval not only improves the signal-to-noise, but also allows for 
the eff ects of line blen ding or overlap and reduces scatter in the 
results dZhu et al.ll2009l) . ideal for the weaker metallic lines. 

Errors were obtained for the cross correlations through esti- 
mating the variance and hence signal-to-noise in each spectral re- 
gion, using the areas identified though spectral clipping to be pure 
continuum. The covariance matrix across 3 adjacent pixels in these 
regions were calculated and the noise reduction factor was com- 
puted. This was combined with the original estimate for the signal- 
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Figure 1. High resolution SOPHIE echelle spectra mounted on the 1.93-m telescope at the OHP taken on JD=2, 454, 106.309. The top panel shows the original 
spectrum, with the red dashed line representing the continuum level. The bottom panel shows the spectrum normalised to unit continuum. 



to-noise, which accounted for the noise correlated across nearby 
pixels in a re-binned echelle spectrograph. Random amounts of 
Gaussian noise centred on the signal-to-noise value were then 
added to the spectra, which were then cross-correlated with the tem- 
plate spectrum once more. The resulting variance in the cross corre- 
lations and line profiles was calculated and the 3<r errors found for 
each pixel, allowing an error bar to be associated with each point. 

3 BALMER LINES 
3.1 Ha Line Profiles 

Fig. [2] shows the Ha line profiles obtained from SOPHIE from 
2007 January 3 rd - February 10 th . All profiles appear P Cygni 
in nature. The profiles show a blueshifted absorption component 
with the edges extending up to around —250 to —300 km s~ with 
the deepest absorption at ~ —50 km s _1 and strength weakening 
as blueshifted velocity increases. The redshifted emission compo- 
nent shows greater variability in strength: only som e profiles show 
a reds hifted emission component as first noted by ICroswell et alj 
dl987h . with the red edge of the emission ext ending from 200 
to 250 km s _1 , similar to those observed by iBastian & Mundtl 
(1985). The presence of a correlation between the variability of 



the strength of the redshifted emission peak and the width of the 
blueshifted absorption trough is initially not entirely clear, although 
iD'Angelo et al.1 (|2000) suggested a negative correlation. Fig. [2] 
shows no evidence a negative correlation is sustained in these data, 
for example on JD=2,454,1 13.317 there is strong redshifted Ha 
emission, combined with broad blueshifted absorption. 

To identify any quantifiable relationship between the variation 
of the absorption and emission components of the observed P Cygni 
profiles, a light curve was constructed for each pixel (every 0.5 km 
s~ ) in the line profile, which were then tested for 5000 different pe- 
riods between 2-40 days, evenly distr ibuted i n freq uency space, us- 
ing the method initially d escribed by[Scargle ( 1982), and revised by 
iHorne & Baliunasldl986|) . The Lomb-Scargle periodogram is based 
on least-squares sine curve fitting. Consequently, it is an ideal tool 
for these data because in addition to robustness to irregular cover- 
age it also enables longer periods to be explored and is commonly 
used to explore up to the approximate total coverage period of the 
data. The significance of the periods were found to false alarm prob- 
ability, (FAP) of less than 0.01 using Fisher's method of random- 

2 A variation in the minimum of the blueshifted absorption may also give 
rise to this variability. However, the superposition of Ha emission compo- 
nents make this difficult to identify. 
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Figure 2. Variations in the Ha profile of FU Orionis obtained from the 
SOPHIE echelle spectrograph mounted on the 1.93-m telescope at the OHP. 
The 21 spectra cover a time frame of 38 nights, with the Julian dates of the 
spectra -2,454,000 displayed on the y-axis. All spectra have a resolution 
of 4 km s — 1 , are barycentrically corrected, and placed in the rest frame of 
FU Orionis. The shaded region depicts the velocities where periodicity was 
found to be significant below the 0.01 FAP level. 
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Figure 3. Contour periodogram showing the significant periods found in 
the Ha line profiles. The 0.01 FAP corresponds to power >7.68. 



ization Jlinnell Nemec & Nemed Il985t) . The periodogram show- 
ing the significant periods found in the Ha line profile are shown 
in Fig. [3] The periodogram highlights the significant periods with 
FAP<0.01, corresponding to a power > 7.68. These were in two 
distinct regions from —122.0 to —105.5 km s _1 and from —32.5 to 
+85.5 km s _1 , shown by the shaded regions in Fig. [2] The range 
of significant periods detected was from 11.59-16.29 days, with 
the most probable value found from the 'centre of power' of all 
the significant powers detected in th e line profile at 13 .28 days. 
This period is similar to that found by Herb ig et alj d2003l) . who de- 
tected a period in the 'fast wind' in the equivalent width, (EW) of 
the P Cygni absorptio n between -110 to -270 km s _1 of 14.847 
days. However, unlike iHerbig et ail d2003l) . we detect the period in 
the line profile itself, and not in the EW of the absorption compo- 
nent. The overlap in vel ocity space of the significant periods de- 
tected with the results o f IHerbig et al] 12003) confirms the stability 
of the period in the blueshifted absorption component of Ha over 
a decade. In addition, Fig.[3]extends the detection of th e significant 
period to less negative velocities than the results of IHerbig et alj 
( 2003). The absence of periodicity from —105.5 to —32.5 km s~\ 
suggests the periodicity detected < — 100 km s _1 is due to vari- 
ations in the FWHM of the blueshif t ed ab sorption component of 
Ha as proposed by iD'Angeloetafl (2000). This 'absent' region 
in Fig. [2] corresponds to the deepest blueshifted absorption on all 
nights observed. Likewise, Fig.|2]also demonstrates less variability 
around the minimum of the absorption trough, w hich could be at- 
tributed to the suggestion bv lCroswell et alj d 1987t) that the Balmer 
lines are so optically thick that they are saturated, and consequently 
are not very vulnerable to variations in physical parameters such as 
density/temperature fluctuations or changes in the mass-loss rates. 
However, at higher negative velocities, further out in the expanding 
wind, the Ha line may no longer be saturated, and environmen- 
tal variations are reflected in the variability of the line profiles in 
such regions. Alternatively, the phase-complexity of the superposi- 
tion of different variable components contributing to the Ha profile 
may prevent any periodicity in this region from being detected. The 
periodicity in the strength of the redshifted emission component is 
the same as the periodicity in the blueshifted absorption compo- 
nent, shown by the alignment in frequency of the significant peaks 
in the periodogram in Fig. [3] This is the first t ime such correla- 
tion h as been detected. Previous observations bv lVittone & Erricol 
have seen periodicity in the redshifted emission component 
of 6.70 days, which is not evident in our results. The dense time 
coverage of our results confirms this period is no longer present, 
and could be a result of aliasing in the earlier observations. There- 
fore, in contrast to previous beliefs, the periodic variability in FU 
Orionis seems to be stable over hundreds of periods, because the pe- 
riod detected here agrees well wit h the significant periods of 13-18 
days found bv lHerbig etal] d2003l) . 

The phase diagram in Fig. [4] shows the phase relationship be- 
tween the periodic behaviour detected in various spectral regions 
phase folded with the variance weighted average period of 13.48 
days from all line profiles investigated. In each case the light curve 
is constructed at the power weighted centre of velocity for the pe- 
riodic region concerned. The 2 nd panel from the top shows the 
two periodic regions of the Ha profile. The amplitude of vari- 
ation corresponding to the oscillations in redshifted emission in 
the Ha profiles (the red line) is greater than those at the faster 
wind velocity (the blue line). This is expected from the appear- 
ance/disappearance of the Ha emission peak, in contrast with the 
broadening-narrowing FWHM of the blueshifted absorption in Fig. 
[U Fig.[4]confirms the 'absent' region in Fig.[2]is not simply lacking 
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Figure 4. Phase diagram comparing the periodicity found in different spectral lines (top five panels), with the periodicity found in the cross correlations, 
(bottom three panels). The spectral lines are phase folded with the variance weighted total average period of all five regions, which is 13.48 days. The spectral 
lines are plotted at the power weighted average velocity of the relevant significant detections (from top to bottom), H/3 at +3.5 km s — 1 , Ha, (in redshifted 
emission at +10.0 km s~ 1 plotted with diamonds, and in blueshifted absorption at —109.5 km s _1 plotted with crosses), Nal D, (with Na Di A5895.92 
at -10.0 km s _1 plotted with crosses and Na D 2 A5889.95 at -11.0 km s _1 , plotted with diamonds), Lil A6707 at -35.0 km s" 1 and Fell A5018 at 
-39.0 km s" 1 . The regions cross correlated are shown in the bottom three panels. The cross correlation profiles were phase folded with a variance weighted 
combination of the individual periods, which averaged 3.64 days. The cross correlations are also plotted at the power weighted average velocity which are, 
(from top to bottom) -9.0 km s" 1 in 6320-6440A, -7.5 km s" 1 in 6138-6191A (from cross correlations with observed template of /3Aqr) and —12.0 km 
s in 5540-5640A (from cross correlations with synthetic stellar spectrum of GO star at T e fj=6000 K, log(g)= 3.0). The solid lines show the sinusoidal line 
of best-fitting for each individual detection, which are plotted in blue when the period is detected at blueshifted velocities (the fit to Na D2 is plotted in light 
blue for clarity), and red corresponding to detections in redshifted velocities. All error bars are shown to the 3 sigma level. 
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Figure 5. Same as Fig.|2]except for the H/3 profile at 4861.34 A. 



of periodicity due to complex superposition of the different variable 
components, since this would require the relative contributions to 
be equal and opposite, which is clearly not the case. The satura- 
tion of the Ha line between —105.5 to —32.5 km s _1 seems a 
more plausible explanation. The sinusoidal curve fitted to the fast 
blueshifted data also appears to have a slight lag of ~ 1.80 days 
in comparison with the sinusoidal curves fitted to the data corre- 
sponding to oscillations in the redshifted emission peak. This could 
stem from the different formation locations of the two line com- 
ponents, with the redshifted emission formed in the disc chromo- 
sphere, and the blueshifted absorption formed further out in the 
expanding wind. 



3.2 H/3 Line Profiles 

The H/3 line profiles obtained from SOPHIE are shown in Fig. 
[5] All profiles show broad asymmetric absorption troughs, with 
stronger absorption on the blueshifted edge, which extends out to 
faster velocities than the redward absorption wing. This is to be 
expected for the H/3 profile, whic h is thought to b e produced in 
an accelerating, massive outflow dHartmanrJ 1 1 9941) . However, all 
the profiles appear more skewed towards blueshifted velocities than 
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Figure 6. Contour periodogram showing the significant periods found in 
the H/3 line profiles. The 0.01 FAP corresponds to power > 7.02. 



previously observed in 1992 bv lHartmann & Calvetl ( fl99l) . which 
may imply an increase in the mass-loss rate since this epoch. The 
maximum depth of the absorption trough appears to remain consis- 
tent at ~ —25 km s _1 throughout our observations. However, there 
is still significant night-to-night variability visible in Fig. [5] For in- 
stance, at slower negative velocities, and low positive velocities the 
gradient of the profile appears highly variable. This could be an 
effect of the varying relative absorption strengths between the red- 
shifted and blueshifted components contributing to the overall line 
profile. In addition, the FWHM of the blueshifted absorption com- 
ponent appears to vary significantly throughout the observations, 
which seems to correlate with the apparent broadening of the Ha 
absorption discussed in section |3~T1 This is evident when compar- 
ing observations on JD=2,454,1 13.317 and JD=2,454, 132.307, dis- 
played in Figs. [2] and [5] both of which show relatively broad and 
narrower blueshifted absorption components respectively. 

The presence of any periodic variability in the H/3 line profiles 
was investigated between 2-40 days with the analysis described in 
section [3~T1 The region of frequency-velocity space identified as sig- 
nificant to below the 0.01 FAP, which represents a power > 7.02 
in these data, corresponds to a period from 8.82 to 14.68 days at 
velocities of —8.0 to +11.0 km s _1 in the line profile. This feature 
can be seen as the red blob in Fig. [6] with a power-weighted aver- 
age period of 12.89 days and at a average velocity of +3.5 km s _1 . 
This is the first detection of periodicity in H/3. It is surprising that 
the periodicity is not seen at faster blueshifted velocities like Ha, 
given the formation of the both lines within the massive outflow. 
However, the examination of the top two panels in Fig. [4] shows 
that the periodicity in the redshifted emission of Ha and redshifted 
side of H/3 appear to be very similar in-phase. This could suggest 
that it is the redshifted emission component in H/3, formed in the 
disc chromosphere, and not a variation in strength of the superim- 
posed H/3 redshifted absorption, which is the source of this peri- 
odicity. On examination of Fig. [5] at the velocities where a period 
was identified (the grey shaded region between —8.0 to +11.0 km 
s" 1 ), it appears this cyclic variability corresponds to the shallow- 
ing/deepening of the 'ridge' connecting the blueshifted absorption 
to the redshifted absorption. Consequently, a variation in H/3 emis- 
sion could produce this effect, and a narrower formation region of 
H/3 emission in comparison to Ha emission within the chromo- 
sphere, or the additional superimposed variation in the redshifted 
absorption component of H/3 may prevent this periodicity being de- 
tected at faster redshifted velocities alike Ha. 
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Figure 7. Same as Fig.[2]except for the Nal D profile of FU Orionis. The 
lower axis is centred around Nal Di at 5895.92 A, and the upper axis is 
plotted with respect to Nal D 2 at 5889.95 A. 



3.3 H7 Line Profiles 

The presence of periodic variability was also tested analogously in 
the H7 profiles. In this case, the observations showed no evidence 
for periodicity in the profiles from —300 to +300 km s" 1 above 
the random noise level, which is higher in this shorter wavelength 
region. 



4 NA I D PROFILES 

The Nal profiles obtained from SOPHIE are shown in Fig. 
[7] All profiles show complex absorption components, with no 
noticeable emission c omponents, as first observed in 1981 by 
iBastian & Mundt; ( fl985l) . The blueshifted profiles both seem to con- 
tain two dominate absorption components, a high-velocity broad 
feature, superimpose d with a narrow, low-velocity absorption dip 
dCroswell et aljfl987l) . The broad blueshifted feature appears to be 
variable during our observations with the edge extending out to 
velocities between —150 to —250 km s _1 , (slightly slower than 
the expansion velocities seen in Ha in Fig. [2j. This is consis- 
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Figure 8. Contour periodogram showing the significant periods found in 
the Na I D line profiles. The 0.01 FAP corresponds to power > 6.67. 



tent with the lines' formation in a differentially expanding wind, 
similar to the broad blueshifted absorption of Ha, but in a more 
compact region, with the wind continuin g to expand throughou t 
the more extende d Ha formation region jBastian & Mundtl l 19851) . 
IBastian & Mundtl d 19851) went on to conclude that the tempera- 
ture of the accelerating region must increase outwards, because 
Ha ge nerally needs higher temperatures to be formed than Na I 
D, yet ICroswell et al.l d 19871) showed the Ha absorption could be 
formed at lower temperatures with an increased mass-loss rate. In 
addition, the relatively broad absorption feature of Na I D seems 
to be comprised of two separate variable absorption components, 
which are visible separately on some nights and appear blended 
together on others, (see JD=2,454,104 319 and JD=2 4 54,13 0.292 
on Fig. [7J. This is cons istent with ICroswell" et al.1 dl987t) and 
Hartmann & Calvetl ill 9951) . who observed a single-component and 



Bastian & Mundt 



19851) who observed two-components in the 
broad absorption features. On the other hand, the narrow low- 
velocity component ~ —3 km s _1 does not appear to vary sig- 
nificantly during our observations, which is consistent with it's pro- 
posed formation in an expanding circ umstellar shell, at large dis- 
tances from the object, as discussed by Bastian & Mundt ( 1985). 

The presence of periodicity in the Na I D line profiles between 
2-40 days were analysed using the same methods described in sec- 
tion |3. 1 1 The significant periods found in this doublet are shown by 
the two red regions in Fig. [8] which are below the 0.01 FAP, with 
power > 6.67. In the Na Di line at 5895.92 A, periodicity is sig- 
nificant to less than the 0.01 FAP between —16.0 to —5.5 km s 
for periods of 12.41 to 17.87 days. The range of significant peri- 
ods detected in the Na Di line not only overlaps with the periods 
detected to less than the 0.01 FAP in the Na D 2 line at 5889.95 A, 
but are also present at the same relative velocities in the line profile, 
shown by the shaded regions in Fig. [7] Significant periodicity was 
detected in Na D2 between —25.0 to —5.5 km s _1 with periodic- 
ity ranging from 11.59 to 17.54 days. This is the first detection of 
periodicity in the Na I D profile and the symmetry of these two indi- 
vidual detections justified their combination to calculate an overall 
power-weighted average period of 14.27 days for the Nal D pro- 
file. There are also some other detections in Fig. [8] that appear to 
have power values above the power corresponding to the 0.01 FAP 
level, however these detections are situated on the very edge of the 
periodogram at small frequencies and do not coincide with an ex- 
tended feature. Therefore these 'significant detections' are assumed 
to be a result of edge affects in the frequency sampling and are not 
considered in further analysis. 
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The origin of the periodicity highlighted by the grey shaded 
regions on the line profiles in Fig. [Vj is difficult to distinguish due 
to the complexity of the blueshifted absorption components con- 
tributing to both lines. However, the periodicity of both lines in 
the doublet appears to be cut off at the edge of the narrow low- 
velocity blueshifted absorption component. This is consistent with 
the formation of this low-velocity component in a slowly expand- 
ing shell at a large distance from the object. Therefore, the peri- 
odicity is more likely associated with the lower velocity end of 
the broader, more variable blueshifted components, which are at- 
tributed to the expanding wind. This is supported by the periodici- 
ties overlapping with the significant detections in Ha and H/3. The 
detection of these periods at relatively low velocities in comparison 
to the maximum velocity of the outflow ~ —150 to —250 km s _1 
in this case, suggests the origin of this periodicity is lower down in 
the expanding wind, in contrast to Ha in section [3~T1 where peri- 
odicity in the blueshifted absorption component was seen centred 
at —109.5 km s _1 . This disparity could be due to a combination 
of the more extended formation region of H a and the relatively re - 
duced temperature sensitivity of the Na I D, dCroswell et ai]|l987h . 
which may prevent a significant periodic signal at the edges of the 
formation region, in the expanding wind. Fig. [4] supports this since 
the 3 rd panel from the top shows the periodicity in both lines phase 
folded with a variance weighted average overall period of 13.48 
days. Both Na I D lines, (Di shown in dark blue, D2 shown in light 
blue) appear to have a slight lag with the Ha emission, thought to 
be formed in the disc chromosphere, and a slight lead over the Ha 
absorption, formed further out in the expanding wind, shown in the 
panel above in Fig. [4] 



5 METALLIC LINE PROFILES 

Line profile variations in Li I A6707 from 2007 January 3 rd to 
February 10 th are shown in Fig. [9] These profiles support the 
accretion disc hypothesis o f FU O ri onis objects as d escribed by 
lHartmann & KenvorJ Jl985t Il987ah ; Izhu et alj (2009), since the 
profile display s a double peaked nature, as initially observed by 
Herbig] dl989h : however variability only exists in the blueshifted 
absorption. Unlike Ha, as described in section 13.11 the width 
of the blueshifted absorption appears to remain constant. The 
blueshifted absorption instead varies in strength, causing the dou- 
ble peaked line profile to appear more and less asymmetrical, see 
JD=2,454,121.333 and JD=2,454,127.377, on Fig.[9] These obser- 
vations coincide with the dates when the Ha profile had no red- 
shifted emission, and strong redshifted emission respectively (see 
Fig-HJl. This could indicate a link between the origin of the variabil- 
ity in the Ha line and Li I A6707. 

The possibility of periodicity between 2-40 days in the Li I 
A6707 line profile was investigated via the method described in sec- 
tion |3. 1 1 producing significant detections at velocities from —43.0 
to —27.0 km s -1 , with periods between 12.21 and 19.78 days, re- 
sulting in a 'centre of power' period at 14.81 days. This is the first 
detection of periodicity in Li I and can be clearly identified as the 
large orange-red structure in Fig. [10] which displays the resultant 
periodogram for the frequency-velocity space analysed. There ap- 
pears to be an addition feature on the periodogram in Fig. 1 101 that 
stands out in yellow above the background at velocities between 
—69.5 to —68.0 km s _1 corresponding to a period between 21.22 
and 29.39 days, and a power-weighted average period of 24.76 
days. This feature is significant below the 0.01 FAP level, which 
corresponds to a power > 6.76. On inspection of Fig. QJJ] it is 
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Figure 9. Same as Fig.f2]except for the Li I profile at 6707.89 A. 
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Figure 10. Contour periodogram showing the significant periods detected 
in the Li I line profile. The 0.01 FAP corresponds to power > 6.76. 



unclear whether this periodicity is connected to the larger feature 
~ —35 km s , as it is located at slightly more negative veloci- 
ties towards shorter frequencies. However, since there is no overlap 
in frequency space between the significant periods detected it is 
treated as a separate detection and for the remainder of this discus- 
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Figure 11. Same as Fig.[2]except for the Fe II profile at 5018.60 A. 
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Figure 12. Contour periodogram showing the significant periods detected 
in the 5018 Fe II line profile. The 0.01 FAP corresponds to power > 6.73. 



sion is classed as noise fluctuations, until further detections prove 
the contrary. 

The location of the main spike in power on Fig.llOlcoincides 
with the deepest absorption in the Li I A6707 profile in Fig. [9] lo- 
cated at ~ —35 km s - , with the connection highlighted by the 
shading of the periodic region. This confirms the periodicity in Li I 



is a result of the varying strength of the blueshifted absorption peak, 
and not a periodic broadening of the line profile. The confinement 
of periodicity to purely negative velocities indicates this could be 
another periodic wind signature, although it was previo usly thought 
the L i I absorption line was not formed in the wind aCalvet et alj 
Il993h . Despite the apparent difference in nature of the periodicities 
on the blue wing with Ha, the periodicity observed is similar in na- 
ture to that of Na I D in Fig. [7] but observed at slightly higher veloc- 
ities. This difference could be down to the separate formation loca- 
tion of the lines in the expanding wind, since the periods are similar. 
The 4 th panel down in Fig.|4]shows the data in the power weighted 
average velocity channel of the significant period detected, (—35.0 
km s _1 ) phase folded with the overall average period of 13.48 days. 
(The range of significant detections of periodicity in Li I includes 
the 13.28 and 14.27 day periods found in Ha and Na I D profiles re- 
spectively, so it is plausible that all these periods could be produced 
via the same mechanism.) The diagram in Fig.|4]shows clear cyclic 
variation and when fitted with a sinusoid, shows the variations are 
roughly in-phase with the variations of the Na I D components, dis- 
played in the panel above. The intermediate phase of both Li I and 
Nal D between that of the chromosphere, (red lines in Fig. |4]> and 
the higher velocity blueshifted components of Ha are consistent 
with both the presence of periodicity and the lines' formation lower 
down in the expanding wind. 

Another metallic absorption lin e not as opticall y thick as Ha is 
the Fe II line at 50 1 8 A, suggested by Herbig] j 19661) to be produced 
further out in the wind in a 'shell' feature of FU Orionis displaced 
~ —80 km s _1 . Fig.ll llshows the A5018 line profiles, which show 
no indication of double peaked structure, but instead are predom- 
inately asymmetric with maximum absorption depth ~ —50 km 
s~ , slightly greater than Li I, more similar to Ha and H/3. On 
some nights, most notably JD=2,454,121.333, the absorption forms 
a sharp dip, which could be attributed to a shell-like sig nature in the 
wind. However, this spectral feature was identified bv lCalvet et alj 
dl993h to be present when the disc rotational velocity is comparable 
to the expansion velocity in the wind, thus the 'shell' signature is 
disingenuous. Consequently, the loss of the region of sharp decline, 
and the production of a more symmetric blueshifted absorption pro- 
file on Fig.[TT]by JD=2,454,127.377, has implications for the mass- 
loss rate and expansion velocity in the wind. The blueshifted ab- 
sorption appears to extend out to ~ —125 km s _1 , however both 
the width and the strength of the absorption profile in general ap- 
pears to vary, comparable to the Ha and Li I A6707 profiles respec- 
tively. 

The periodogram shown in Fig. [T2] was obtained via the 
method described in section [3~T1 over the same parameter space 
for comparison. It shows less significant detections of a period than 
the Li I A6707 profile, but there are still periods significant to the 
0.01 FAP with velocities from —43.5 to —24.0 km s _1 with values 
between 12.31 and 15.68 days. This corresponds to the largest fea- 
ture in Fig.[T2] which has a power weighted average period of 13.57 
days. The dominant periodic feature detected in this line profile is 
yet again of a similar value to that of Ha, H/3 Na I and Li I, pro- 
viding further evidence that there is a wind period in FU Orionis. 
This is confirmed in Fig. [4] where the 4 th panel from the bottom 
shows the Fell A5018 data in the power weighted average veloc- 
ity of the periodic feature (—39.0 km s _1 ), phase folded with the 
average overall wind period of 13.48 days. The cyclic variations 
are fitted with a sinusoid, (solid blue line), and the fit is better than 
that for Na I D and Li I in the two panels above. However, unlike 
the previous line profiles, the phase of the Fell line does not lie 
within the range expected for formation in the expanding wind, and 
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Figure 13. Reference spectra used in the cross correlations with FU Orionis, shown in the 5540-5640 A region. The bottom black line is a synthetic stellar 
spectrum of T c g = 6000 K, log(g)= 3.0, the blue line is the same spectrum smoothed with a Gaussian with FWHM of 0.37 A, which was found to 
correspond to the FWHM of the unblended, unsaturated lines in the j3 Aqr spectrum. The top red line is the spectrum of Aqr, obtained from the SOPHIE 
echelle spectrograph mounted on the 1.93-m telescope at the OHP with a resolution of 4 km s _1 . 



seems to respond even before the chromospheric lines in the top 
two panels in Fig. [4] Although there are less significant detections 
of a period to less than 0.01 FAP visible in Fig. 1 121 than in Li I 
A6707 (see Fig. 1 10b. p ossibly due to the high temperature sensitiv- 
ity of the Li I line dCalvet et al.fl993l) . the amplitude of variation is 
greater in the Fell line at A5018 (see Fig. |4j- On examination of 
Fig.Qjj the velocities at which this period was detected, (shown by 
the shaded region) seem to correlate with the appearance and dis- 
appearance of the absorption 'dip', confirming the assumption of 
the period's origination in a wind. The fact that there are less sig- 
nificant detections of periodicity could indicate that the effects of 
the mechanism for producing this period in the wind may become 
reduced at greater distances from the disc where the Fell A5018 
absorption line is thought to originate. 



As in the Li I A6707 profile, there is also another feature vis- 
ible on the periodogram in Fig.[l2]at slightly lower frequency and 
more negative velocity than the dominant feature, corresponding to 
a weak detection of periodicity, above the 0.01 FAP level, (power 
> 6.73). The range in velocity of this yellow feature is —128.0 
to —123.5 km s~ and spans periodicities between 19.52 to 23.91 
days. This is comparable to the detection discussed in Li I A6707, 
where the authenticity of the detection was unclear. In Fe II A5018 
this detection has a power weighted period of 21.66 days, slightly 
lower than the comparable signal in Li I A6707, yet there is a sub- 
stantial overlap in the significant periods detected in both lines. The 
velocities at which the period manifests itself are slightly higher in 
Fell A5018 than in Li I A6707, however this could be ascribed to 
the different in formation locations in an expanding wind. The re- 
currence of this feature at a comparable periodicity to an ambigu- 
ous detection in Li I A6707 that affects a larger area of surrounding 
parameter space in Fig.[T2]then inQJJ] makes this detection slightly 
less likely due to noisy line profiles and may in fact be linked to the 
main source of periodicity in both cases. 



6 CROSS CORRELATIONS 

Some of the metal lines were too weak and too noisy to anal- 
yse the line profile alone. In order to increase the signal-to-noise 
ratio, these profiles were cross correlated with a template spec- 
trum of P Aqr and a synthetic stellar spectrum of a GO star with 
T e fi = 6000 K, log(g)= 3.0 to identify any potential effects of rel- 
ative line shifts, as described in section [2] The standard spectra in 
the 5540-5640 A spectral region are shown in Fig. [T3] where the 
bottom black line is the synthetic spectrum and the blue line shows 
the Gaussian smoothed synthetic spectrum to match the FWHM 
of the unblended, unsaturated lines in the observed spectrum of /3 
Aqr, shown by the top red line. The regions cross cor related were 
chosen in order to compare with previous results by [Herbig et alj 
where a 3.542 day period ha s already been identified in th e 
centroid of the cross correlation, and Hartmann & Kenvo^ <1987iJ) . 
who found variability in their cross c orrelation profile of the 6170 
A region. The regions investigated bv lHerbigetai] ( l2003l) were the 
6320-6440 A region and the 5540-5640 A region, which were ini- 
tially identified in order to avoid lines with wind components. The 
normalised SOPHIE echelle spectra taken from 2007 January 3 rd to 
February 10 th are shown in Figs.[l4][T5]and[T6] The spectra in the 
5540-5640 A region, (see Fig. [14} were normalised using a lower 
bias level (0.97 instead of 0.98), because the continuum appeared 
under-subtracted when the higher value was used. All spectra show 
evidence of weak absorption lines and some weak lines appear dou- 
bled. This is most notable in Fig.ll5lwhere the 6162.17 A Cal and 
6169.56 A Cal absorptions appear doubled on some nights, indi- 
cating a disc orig in, although the pro fi les are highly variable. Other 
lines noted by lHartmann & KenvorJ J1987al) to show evidence of 
double peaked structure in this spectral region are the 6147.8 A Fe I 
and Fell line blend, 6151.62 A Fel, 6157.73 A Fel, and 6166.44 
A Cal. There does appear to be some evidence of a skewed ab- 
sorption profile in the 6141.73 A Fel line, indicating the presence 
of mass-loss. Cross correlating these spectral regions provides a 
distinct advantage as not only does it produce an average profile 



© 2012 RAS, MNRAS 000,mi20l 



12 S. L. Powell et al, 

142.333 



X! 

w 
E 



CNJ 




137.335 
132.307 
131.304 
130.131 
129.307 
128.307 
127.377 
121.333 
120.317 
1 19.420 
1 18.320 
1 16.372 
15.310 
1 15.287 
1 14.316 
1 13.317 
109.317 
107.320 
107.301 
106.309 
105.334 
104.319 



5540 5560 5580 5600 5620 5640 

Wavelength (A) 

Figure 14. Spectra of FU Orionis in the 5540-5640 A region obtained from the SOPHIE echelle spectrograph mounted on the 1.93-m telescope at the OHP. 
The 21 spectra have a resolution of 4 km s~ 1 and cover a time frame of 38 nights, with the Julian dates of the spectra -2,454,000 displayed on the right hand 
y-axis. All spectra have been barycentrically corrected, normalised to unity, placed in the rest frame of FU Orionis, continuum subtracted with a bias level of 
0.97 and shifted to show night to night variability. 
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of 0.98 for the continuum subtraction. 
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Figure 16. Same as Fig. [14] but for the 6320-6440 A region using ; 
bias level of 0.98 for the continuum subtraction. 
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weighted according to the standard line strengths, but it also allows 
for lines which overlap and blending, abundant in these regions, 
which is clearly a problem for identifying the weak lines in Figs. 
[14] [15] and [T6] The cross correlation functions (CCFs) of each night 
were tested for periodicity at 5000 different frequencies between 2- 
40 days to quantitatively compare the results with the line profiles 
in sections[3][4]and[5] 

The CCF of the 6138-6194 A region, displayed in Fig. [HI 
showed the clearest detection of a period significant to less than 
0.01 FAP, (power > 6.88) in velocity channels -14.0 to -1.5 km 
s _1 . The significant periods range from 3.52 to 3.70 days and have 
a power weighted average of 3.61 days. This region is also clearly 
visible above the noise in Fig. [T8] This period is similar in value 
to that detected previously, yet it is distinctly different in signature, 
corresponding to a cyclic strengthening/weakening of th e blue wing 
of the line profiles in this region. To compare with iHerbig et alj 
( 2003), the radial velocity of centroid of the CCF was found for 
each observation and tested for periodicity between 2 and 40 days, 



yet no significant period was found below the 0.01 FAP. This was 
unexpected, as the periodicity has only been observed on the blue 
wing of the CCF. To make a more direct comparison , the spectra 
were d egraded to a resolution of 13 km s _1 used bv lHerbig et alj 
(2003). At this resolution a longer period of 13.12 days, (compa- 
rable to the periodicity previously observed in the line profiles di- 
rectly) was detected in the centroid, (or weighted average velocity) 
of the cross correlation in the 6138-6194 A region, which was sig- 
nificant to the 0.01 FAP level. This could possibly indicate the ef- 
fects of the wind are more significant in this spe ctral region, than in 
the regions investigated by Herb ig et al.1 d2003l) . and the shorter pe- 
riodic effects may be masked in the centroids due to the other vari- 
able features in the CCF, most notably the effects of the wind. The 
power weighted average radial velocity of the 'blob' in Fig.ll8lwas 
— 7.5 km s -1 , and the phase folded data at this velocity is shown 
in the 2 nd panel from the bottom in Fig. [4] The precise velocity at 
which this periodicity is detected is difficult to validate since spec- 
tral lines in the observed GO lb template of /3 Aqr may be shifted 
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Figure 17. Cross correlations of FU Orionis in the 6138-6194 A region with a reference spectrum of f3 Aqr obtained from the SOPHIE echelle spectrograph 
mounted on the 1.93-m telescope at the OHP. The cross correlation function is the black line, with the upper and lower 3cr errors shown by the dark blue and 
light blue lines respectively. When two observations were taken on a single night a dark green line is used to represent the second observation. The average 
over all cross correlations in this region is shown by the red dashed line. 
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Figure 18. Contour periodogram showing the significant periods found in the 6138-6194 A region. The 0.01 FAP corresponds to power > 6.88. 
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Figure 19. Same as Fig.rfT] but for the 6320-6440 A window. 
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Figure 20. Contour periodogram showing the significant periods found in the 6320-6440 A region. The 0.01 FAP corresponds to power > 6.73 
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Figure 21. Same as Fig. 1 171 but for the 5540-5640 A window, which was continuum subtracted with a lower bias of 0.97 and cross correlated with a GO 
synthetic stellar template at T c g = 6000 K, log(g)= 3.0. 
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Figure 22. Contour periodogram of the 5540-5640 A cross correlation, 
which was continuum subtracted using a lower bias of 0.97 and cross 
correlated with observed reference template of j3 Aqr. There are no de- 
tections in the cross correlation function below the 0.01 FAP which cor- 
responds to power > 5.88. 
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Figure 23. Contour periodogram showing the significant periods found in 
the 5540-5640 A cross correlation, which was continuum subtracted with 
a lower bias level of 0.97 and cross correlated with a GO synthetic stellar 
template at T cff = 6000 K, log(g)= 3.0. The 0.01 FAP corresponds to 
power > 6.16. 
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due to nightly fluctuations in air mass, which could ultimately shift 
the velocity peaks in the cross correlations. However, the veloc- 
ity range of the periodicity was reinforced by cross correlations 
with the synthetic stellar spectrum at T e ff = 6000 K log(g)= 3.0, 
which confirmed periodicity between 3.54 to 3.71 days to below 
0.01 FAP spanning velocities from —16.0 to —4.0 km s -1 , with a 
power weighted average period of 3.62 days. 

The detections below the 0.01 FAP level, (power > 6.73) in 
the 6320-6440 A region cross correlated with the observed tem- 
plate of /3 Aqr were periods from 3.56 to 3.71 days, spanning 
velocity channels from —18.5 to 0.0 km s _1 , which produced a 
weighted average period of 3.63 days. Fig. [T9] shows the result- 
ing CCFs, which look very similar to the cross correlations of the 
6138-6194 A region in Fig. [17] Fig.l20lshows the significant detec- 
tions with the red blob at the centre of the figure, at a frequency 
of 0.28 days -1 . When the centroids of the CCFs in this region 
were tested for periodicity, yet again there was none found to be 
significant to below the 0.01 FAP level. However, when the spec- 
tra and observed te mplate were degraded to the resolution used by 
Herbi get alj d2003l) in this region, the centroids of the cross corre- 
lation peaks were found to have a periodicity of 13.69 days signifi- 
cant to the 0.01 FAP level. This could possibly indicate an increase 
in the mass-loss rate, which would manifest wind signitures more 
readily in this epoc h. This may also be a possible explanation why 
lHerbigetalJd2003l) noted a poor fit to the sinusoid for the observa- 
tions taken in the year 2000. The velocities at which this periodic 
signal was detected were confirmed by the cross correlations with 
the synthetic stellar spectrum, which also detected a period of 3.63 
days between —19.0 to —1.0 km s _1 , significant to below the 0.01 
FAP. The phase folded data at —9.0 km s -1 fitted with a sinusoid 
shown in the blue solid line, can be seen on the 3 rd panel from the 
bottom in Fig.|4]with error bars plotted to the 3 sigma level. 

Our analysis showed no significant periodicity below the 0.01 
FAP corresponding to a power of 5.88 in the 5540-5640 A region 
when cross-correlated with f3 Aqr. Unsurprisingly, no period was 
detected in the centroid of the cross correlation of the spectra and 
observ ed template degraded to the same resolution as lHerbig et alj 
d2003l) in this spectral region. However, a period between 3.63- 
3.65 days was found in the same spectral region when cross corre- 
lated with the synthetic template, which was significant below the 
0.01 FAP, corresponding to powers above 6.16 in this case. The 
weighted average period was found to be 3.64 days, from —15.0 to 
—9.0 km s -1 and centred at —12.0 km s~ . The resulting cross cor- 
relations with the smoothed synthetic spectrum are shown in Fig. 1211 
with the periodograms displaying the null detection using the /? Aqr 
template, and the significant periods found with the synthetic tem- 
plate in Figs. [22] and [23] respectively. The background level at the 
edges of the CCFs are slightly higher in Fig. |2T] compared to Figs. 
[T9l and [T7J because of the lower bias level used in the continuum 
subtraction. However the shape remains unaffected in this adjust- 
ment, and on visual examination the variability seems to correlate 
with the cross correlations of the 6320-6440 A and 6138-6194 A 
spectral regions using the observed template of /3 Aqr. The absence 
of periodicity in the cross correlations of the 5540-5640 A region 
with the observed template and yet the significant detection with 
the synthetic template, could be manifested by the slightly differ- 
ent relations between spectral line strengths in these two spectra, 
shown with the top two lines plotted in Fig. [T3] Figs. [22] and [23] 
are extremely similar and both show a clear feature around the pe- 
riodicity detected to be significant to below a FAP of 0.01 in the 
cross correlations with the synthetic template. Therefore although 
this feature is not detected to be significant in Fig. [22] to such a 



high level when the observed template of /3 Aqr is used, it is still a 
clear feature above the noise in the periodogram and hence the lack 
of high significance of the feature in this case may be due to the 
slight differences in the line weightings producing slightly differ- 
ent weights in the cross correlation function and consequently, the 
periodogram. In the case of the synthetic standard, these slight dif- 
ferences have contributed to the significance of the cyclic variabil- 
ity in this spectral region. The absence of any significance to less 
than 0.01 FAP with the /3 Aqr standard may therefore not indicate 
the lack of cyclic variability, but rather inadequate signal strength 
in this spectral region possibly due to a different mix/ratio of line 
species in the observed template or a slightly different origin of 
these absorption lines within the disc itself in comparison to the 
spectral regions previously discussed. The bottom panel of Fig. [4] 
shows the data at —12.0 km s _1 phase folded with a period of 3.64 
days and fitted with a sinusoid. For the reasons discussed above, the 
amplitude of variation is small in comparison with the cross corre- 
lations in the other spectral regions, shown in the 2 panels above in 
Fig. g] 

The in-phase nature of the periods found in the cross corre- 
lation functions, along with the extensive overlap in velocity and 
range of significant periods identified to less than 0.01 FAP is 
strong evidence that all these cyclic variations are produced from a 
single origin. Consequently, the three periods identified above were 
combined to find a variance weighted total average period of 3.64 
days, which was used to phase fold the data in all regions plotted in 
Fig. [4] The periodicity could be attributed to a slow moving wind, 
as it is predominantly present at small negative velocities. However 
it does not appear to be related to the wind periods found in the line 
profiles in sections [3] [4] and [5] as it is roughly one quarter of the 
value. 

To complete the analysis of the cross correlations, the FU On- 
onis spectra from each night were cross correlated with the average 
spectrum and the resulting CCFs were tested for periodicity via 
the same analysis. Interestingly in the 6138-6194 A region, this re- 
sulted in periods between 12.31 to 14.32 days, significant to below 
0.01 FAP with a power weighted average of 13.19 days, which are 
similar to the periodicities de tected in the centrio ds of the spectra 
degraded to the resolution of iHerbig et all d2003h and in sections 
[3][4]and[5] Additionally, this periodicity was detected at velocities 
between —80 to —24 km s , more comparable to the fast veloc- 
ities of an expanding wind. No such periodicity was detected in 
the other spectral regions, p ossibly due to a dif ferent mix of line 
species and their selection by Herbig et al. ( 2003) to avoid line with 
wind components. There was no evidence of the 3.64 day period 
in this analysis. However, the resulting cross correlations are ex- 
tremely broad due to the large FWHM of spectral lines in the FU 
Orionis spectra in comparison with the FWHM of the standard tem- 
plates, (see Figs. 1 131 and 1 141 for comparison). Therefore the CCFs 
produced using the average spectra as a template are insensitive to 
small changes around the line centre, which would be necessary to 
detect the 3.64 day period. Consequently, there appears to be more 
than one mechanism producing asymmetries in FU Orionis. 



7 DISCUSSION 

A summary of the velocity channels where periodicity was detected 
to be significant below the 0.01 FAP level in each spectral region 
is shown in Table[2] Our high resolution spectral monitoring of FU 
Orionis has yielded two principle results: 

i) We have confirmed the periodic variability of blue shifted 
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Table 2. Summary of periodicities detected - ( 1 ) Line/spectral region where 
periodicities were detected; (2) Weighted average velocity of periodic fea- 
ture with the limits showing the full range of velocity channels where peri- 
odicity was significant to below the 0.01 FAP; (3) Weighted average peri- 
odicity detected with the full range of periodicities detected above the 0.01 
FAP shown in the limits. 
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function using an observed template of a GOIb star namely f3 Aqr. b The 
periodicity in this region was detected in the cross correlation function 
using a synthetic spectrum of a GO star with T e g = 6000 K, log(g)= 3.0 
as a reference template. 



Ha absorption claimed by iHerbig et al. I fe003l) . The 'centre of 
power' for periodic variation is at a period of 13.48 days, (which 
we will henceforth term 'the wind period') but the range of periods 
suggested by our anal ysis (see T able [2}, also includes the 14.847 
day period found by Herbig et al.l d2003h . In addition, we have for 
the first time, found periodic variation (at the same period) in the 
redshifted emission component of Ha, the redshifted absorption 
component of H/3, and the blueshifted absorption components of 
Na I D, Li I A6707 and Fell A5018. Although the 'centre of power' 
for periodic variations in these lines are slightly different, (13.27, 
12.89, 14.27, 14.81 and 13.57 days respectively) the range of al- 
lowed periods at a false alarm probability of 0.01 all overlap and 
we will assume that the same wind period, (a variance weighted 
average of 13.48 days) applies to all these lines. 

ii) We have also confirmed the periodic variation of the cross- 
correlation function constructed from spectral regions that, lacking 
strong wind features, are believed to originate from the photosphere 
of the disc. The detected period (around 3.64 days in all the spectral 
windows investigat ed) is close to the value of 3.54 days found by 
Herbi g~et alj d2003l) . We find that variability is largely confined to 
low velocities on the blue wing of the cross-correlation function 
(> —18.5 km s _1 and with a velocity centroid at ~ —9.0 km 
s -1 ); a change i n the characteristics of the outflow in FU Orionis 
may explain why H erbig et al.l d2003l) . with a velocity resolution of 
~ 13 km s _1 compared with 4 km s -1 in the present study, also 
found significant variation on the red wing. 

In what follows we will assume that the 13.48 period is associ- 
ated with modulation of the wind. The fact that the 3.64 day period 
is only detectable on the blue wing of the cross correlation function 
raises the question of whether it might not also be a wind period, 
even though the spectral regions from which the cross correlation 



was calculated are believed to originate in the disc photosphere. We 
are dissuaded from invoking a wind component in these lines on the 
grounds that in any case the period is quite different from the 'wind 
period' described in i) above. We note that the two periods are con- 
sistent with being in the ratio of 1 : 4 although the range of allowed 
periods at a false alarm probability of 0.01 or below does not allow 
us to assert this with any confidence. Consequently, we will assume 
that there is no necessary connection between the phenomena giv- 
ing rise to these two periods and that they are associated with two 
physically distinct regions, the wind and the disc photosphere. We 
will briefly revisit this assumption at the end of this section. 



7.1 The periodicity of the wind 

We have found clear evidence (see the upper four panels of Fig. [4) 
that the variations in the redshifted emission of the Ha and H/3 lines 
are in-phase with each other, with the blue shifted absorption of the 
Na I D and Li I A6707 displaying an intermediate phase. The vari- 
ation of the blueshifted absorption of Ha is however phase lagged 
by about 1.8 days. 

This pattern is consistent with a picture in which the wind is 
modulated at its base in velocity and/or density/temperature. The 
first lines to respond are those formed at relatively low velocities 
close to the wind base while the blue shifted Ha absorption (which 
is at considerably higher velocities) responds after a time equal to 
the flow time across the region of acceleration. This hypothesis al- 
lows us to estimate the scale of the region where the wind is ac- 
celerated up to ~ 100 km s _1 as being the product of 100 km s 
_1 and 1.8 days, i.e. ~ 10 12 cm. It should be noted however, that 
the phase of the Fe II line does not fit this picture. The lack of sig- 
nificant periodic behaviour at still higher velocities could either be 
due to the shearing out of high density pulses in the wind beyond 
this region or else because the higher velocity absorption originates 
from a larger region in which a number of consecutive pulses are 
always simultaneously present. 



7.2 The periodicity of photospheric absorption 

A simple interpretation of periodic variability in the line profiles 
from a disc is in terms of an orbiting hotspot which enhances the 
flux (and hence contribution to the absorption line spectrum) at a 
particular velocity. The problem with this simplest model is that 
it would predict that enhanced absorption should switch between 
the blue and red wings of the line as the hotspot alternatively ap- 
proaches and recedes from the observer. However, we find little 
variability on the red wing of the line and so we need to argue that 
the receding hotspot is in some way concealed from the observer. 
This necessarily involves the invoking of three dimensional effects 
in the disc such as a warp. If for some reason the disc surface is 
distorted in the region of enhanced dissipation then the radiation 
may escape preferentially in the forward direction of the orbiting 
feature. 

So far we have been deliberately unspecific about the nature 
of the hotspot or warp. We bear in mind that our confirmation of 
Herbig's period implies that this is a phenomenon that is stable 
over of order a thousand orbits or more and is thus inconsistent 
with a transient heating ev ent which would be sheare d out over 
many fewer orbital periods. IClarke & Armitagel (2003) suggested 
that a long-lived hotspot would be produced if the inner disc of 
FU Orionis contained an embedded (proto-)hot lupiter, as required 
by the model for the triggering of rapid-rise objects (such as FU 
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Orionis and V105 7 Cygni) suggested bv lClarke&Sverl l ll996|) and 
lLodato & Clarke] (2004). In this case, the 'hotspot' is simply the 
emergence of the energy liberated by accretion on to the planet 
as it diffuses out through the overlying disc. If, in addition, the 
disc is warped in the vicinity of the planet, (as would result if 
the planet and disc planes are not perfectly aligned) then this ra- 
diation will preferentially emerge through the region of the disc 
surface which offers the lowest optical depth to escaping photons; 
thus there are geometries for which the hotspot signature would be 
more detectable when the planet is approaching the observer than 
on the opposite phase. 

In a similar vein, one could instead invoke a magnetised 
hotspot on t he di sc surface, (supported by the detection by 
Donati ct al. (2005)) and again, if the axis of the magnetic field 
is misalign ed with that of the disc, one expects a warped structure 
in the disc JEaHl 1 999l [20031 ; iFoucart & Laill20 1 ll) . The stability of 
the observed period implies that the magnetic topology needs to 
remain constant on time-scales of ~ 10 years. 

Note that it is not a necessary ingredient of the model that the 
hotspot co-rotates with the local disc flow, although obviously in 
the planet case, both hotspot and disc orbit at the local Keplerian 
velocity. If this is the case then a period of 3.6 days corresponds to 
an orbital radius of 7 x 10 11 A/ 1 ^ 3 cm where M is the mass of the 
central star in solar masses. The maximum velocity shift would 
in this case be 140M 1/3 sini km s -1 > 60M 1/3 km s -1 where in 
the latter inequality we have adopted the mi nimum value of sin i 
quoted in the literature, dKenvon et alj|l988h and M is the stellar 
mass in solar masses. Even in the case of rather a low mass star 
M ~ 0.3, the resulting peak projected velocity is around 40 km s 
_1 . Inspection of Figs. [T7lll9| andl21lshow that this is close to the 
maximum velocity at which there is detectable variation on the blue 
wing of the cross correlation function, although a significant period 
is only detectable for velocities less than 18.5 km s _1 . The ampli- 
tude of the velocity variation is thus marginally consistent with a 
Keplerian hotspot but only if one adopts rather low values for the 
stellar mass and inclination of FU Orionis. 

On the other hand, we have more freedom in the case of a mag- 
netically generated hotspot since in this case the concentration of 
magnetic field lines associated with the hotspot are not necessarily 
co-orbital with the disc - it is now possible for the combination of 
low velocity amplitude and short period to be associated with an 
interaction between the field and the disc occurring well outside 
co-rotation. Note that if one identified the rotation period of the 
magnetic field with that of the star, then it would be rotating at a 
velocity modestly less than break-up velocity. 

7.3 The relationship between the periods 

We finally return to the issue of whether there can be any con- 
nection between the two varieties of periodic behaviour described 
above. The phase relationship between the low and high velocity 
wind components suggests that we are witnessing the effects of 
pulsed input at the base of the wind flow rather than the effects of ro- 
tational modulation of an intrinsically steady but non-axisymmetric 
wind. In this case we need some process to modulate the injection 
of momentum into the wind close to its base where it interacts with 
the disc. This process does not share the same period as the modu- 
lation of the photospheric lines. Nevertheless it is possible in prin- 

3 Note that this radius is factor 2 — 3 larger than the b est-fitting inner disc 
radius in FU Orionis employed bv lKenvon et alj Jl988h 



ciple for the warped structure that we have invoked to explain the 
modulation of the photospheric lines may also have some influence 
on modulating the mass uptake in the wind (note the similarity in 
spatial scales between the inferred radius of a Keplerian hotspot 
and the vertical height in the wind from which we infer that the vari- 
able wind signatures originate.) This would require that the warped 
structure, in addition to rigid rotation at the orbital period of the 
putative hotspot, also underwent low amplitude tilt motions normal 
to the disc plane. Clearly we have no detailed calculations to sup- 
port this suggestion but merely point out that this could represent 
an economical solution to the 'two period' problem. 
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